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Abstract 

The Euclidean-time hadronic current correlation functions, Gp{t,T) and Gy(r, T), of pseu- 
doscalar and vector currents have recently been calculated in lattice simulations of QCD and have 
been used to obtain the corresponding spectral functions. We have used the Nambu-Jona-Lasinio 
(NJL) model to calculate such spectral functions, as well as the Euclidean-time correlators, and 
have made a comparison to the lattice results for the correlators. We find evidence for the type 
of temperature dependence of the NJL coupling parameters that we have used in previous stud- 
ies of the mesonic confinement-deconfinement transition. We also see that the spectral functions 
obtained when using the maximum-entropy-method (MEM) and the lattice data differ from the 
spectral functions that we calculate in our chiral model. However, our results for the Euclidean- 
time correlators are in general agreement with the lattice results, with better agreement when 
our temperature-dependent coupling parameters are used than when temperature-independent pa- 
rameters are used for the NJL model. We also discuss some additional evidence for the utility of 
temperature-dependent coupling parameters for the NJL model. For example, if the constituent 
quark mass at T = is 352 MeV in the chiral limit, the transition temperature is Tc = 208 MeV 
for the NJL model with a standard momentum cutoff parameter. (If a Gaussian momentum cutoff 
is used, we find Tc = 225 MeV in the chiral limit, with m = 368 MeV at T = 0.) The introduction 
of a weak temperature dependence for the coupling constant will move the value of Tc into the 
range 150-170 MeV, which is more in accord with what is found in lattice simulations of QCD with 
dynamical quarks. 

PACS numbers: 12.39.Fe, 12.38.Aw, 14.65.Bt 
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I. INTRODUCTION 



Euclidean-time hadronic current correlation functions contain information about the 
hadronic spectrum, including information about the temperature dependence of hadronic 
masses, the widths of these states and the eventual disappearance of such status from the 
spectrum at sufficiently high temperature. We have been interested in recent calculations 
of spectral functions of hadronic current correlators [1-3]. These calculations make use of 
lattice results for the Euclidean-time correlation functions, Gp{t,T) and Gv{t,T), of pseu- 
doscalar and vector currents. For example, for the pseudoscalar correlator, one considers 
the equation 



and attempts to obtain values of the spectral function ap{uj,T) by using the maximum- 
entropy-mcthod (MEM) [4, 5]. (A similar procedure is used to obtain (Tv{uj,T).) The 
MEM method is used, since the inversion of Eq. (1.1) to obtain ap{uj, T) is an "ill-conditioned 
problem" [ 2 ] , given the limitations of the lattice data. 

The necessity of using such advanced numerical procedures is discussed in Ref . [ 2 ] . As 
stated there, one has to introduce assumptions as to the characteristics of the desired solu- 
tions when using the MEM analysis. At this point, the physical significance of the results for 
the spectral functions obtained with the MEM procedure is uncertain in the case of resonant 
structures with energies of several GeV found using that method. 

Recently we have made some calculations of hadronic current correlators and their asso- 
ciated spectral functions [6]. Once we calculate the spectral functions using the Nambu- 
Jona-Lasinio (NJL) model, we may obtain Gp{t,T) and Gv{t,T) by using Eqs. (1.1) and 
(1.2). We may then compare our results for these Euclidean-time correlators with those 
calculated using the data obtained in the lattice simulations of QCD. 

Since we have described our method of calculation of the spectral functions in an earlier 
work [ 6 ] , we here relegate that material to the Appendix for ease of reference. A novel feature 
of our work is use of temperature-dependent coupling parameters for the NJL model. We 




(1.1) 



with 



K{t,u;,T) 



cosh[a;(r- l/2r)] 
sinh(a;/2T) 



(1.2) 
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have made use of such parameters in earher studies of meson properties at finite temperature 
and have described the mesonic confinement-deconfinement transition [ 7 ] . While the use of 
temperature-dependent couphng constants for the NJL model is novel, we may note that the 
QCD coupling constant has been made temperature-dependent in a study of the quark 
-gluon plasma equation of state in Ref. [8] and [9]. (See Fig. 1 of Ref. [9].) Values of the 
pressure, energy density and baryon density were presented as a function of temperature 
and quite good fits to QCD lattice results were obtained in Refs. [8] and [9]. 
In Ref. [ 6 ] we have introduced the spectral functions, 

ap{LU, T) = -ImCp(t,;, T) , (1.3) 

TT 

and 

av{uj, T) = -ImCy (u;, T) . (1.4) 

TT 

The relation between these functions and the ones that appear in the literature, ap{uj,T) 
and ay (cu, T) [ 1-3 ], is given in the Appendix. [SceEqs. (A30) and (A31). ] In the Appendix 
we describe the procedures used in the calculation of lmCp{uj,T) and lm.Cv{oJ,T). 

The organization of our work is as follow. In Section II we present some of the results 
obtained for the functions lmCp{PQ,T) / and lm.Cv{PQ,T) / P^. (We identify Pq with u, 
which is used in the literature for the corresponding quantity.) In Section III and IV we 
compare the values of Gp(r, T) and Gyij, T) obtained in our calculations and in the lattice 
simulations [2]. In those sections we discuss the evidence for the temperature-dependent 
coupling parameters we have used in other works [6, 7]. Finally, Section IV contains some 
additional discussion and conclusions. 

II. EUCLIDEAN-TIME CORRELATORS OF PSEUDOSCALAR HADRONIC 

CURRENTS 

The formalism reviewed in the Appendix is used to obtain the results presented in this 
section. In Fig.l we show the values of ImCp(Po, T)/Pq that were first presented in Ref. [6]. 
The large peak at T/Tc = 1.2 has its origin in the properties of a pion-like mode that 
is present after deconfinement has taken place. We note that no resonance is seen for 
T/Tc — 3.0 [dot-dashed line], but resonance behavior is present for T/Tc — 1.5 [dashed 
line] and T/T^ = 2.0 [dotted line]. 
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FIG. 1: Values of ImCp{Po,T)/P^, obtained in Ref. [6], are shown for values of T/T^ = 1.2 
[solid line], T/T^ = 1.5 [dashed line], T/Tc = 2.0 [dotted line], T/Tc = 3.0 [dot-dashed line], 
T/Tc = 4.0 [double dot-dashed line], and T/Tc = 5.88 [short dashed line]. Here we use Gp{T) = 
Gp[l- 0.17 (T/Tc)] with T^ = 0.150 GeV and Gp = 13.49 GeV-^. 

Before we proceed, it is useful to discuss the properties of the integral in Eq. (1.1). We 
wish to show that the calculation of Gp{t,T) is sensitive to the properties of ap{ijj,T) for 
relatively small u, when tT ~ 0.5. On the other hand, when tT is small (or near 1), the 
integral is sensitive to the values of (7p(t, T) for large a;, where a"p(T, T) increases as o;^ and 
is largely model independent. These features may be seen in Fig.2, where the solid line 
represents our calculated values of ap{T,T) = (l/7r)ImCp(cj, T) for T = 1.5 Tc. The dashed 
line shows K{u!,t,T) for tT = 0.5, while the dotted and dot-dashed lines show K{u!,t,T) 
for tT — 0.05 and tT — 0.10, respectively. We note that the 1/u; singularity in K{uj,t,T) 
for small uj is compensated by the behavior of the spectral function at small a;, (Tp(a;, T) ~ a;. 
These remarks also pertain to the calculation of G'y(r, T) made using cry (a;, T). 

In Fig. 3 we show the data we have taken from Ref. [ 2 ] as squares. The sohd hue shows the 
result of our calculation of Gp(t, T)/T^, as obtained from our calculated values of ap{ui, T) 
aX T — 1.5 Tc. The deviation of our calculation from the data for tT ~ or tT ~ 1 
is inconsequential, given our remarks made in the previous paragraph. We also show the 
result for the case Gp{T) — Gp — 13.49 GeV~^ as a dotted curve. When comparing the 
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FIG. 2: Values of K{u,t,T) are shown for T = 1.5 and tT = 0.05 [dotted line], tT = 0.10 

[dot-dashed line], and tT = 0.50 [dashed line]. The solid line represents (l/7r)ImCp(a;, T). Here, 
we use the notation to = Pq and have put Tc = 0.150 GeV. 
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FIG. 3: Values of Gp(r, T) /T^ are shown as a function of tT, with T = 1.5 T^. The solid line repre- 
sents the result of our calculation made for Gp{T) = Gp [1 — 0.17 (T/Tc)] with Gp = 13.49 GeV~^. 
The dotted line is obtained when we use a constant value Gp{T) = Gp. The data (squares) are 
taken from Ref. [2] for the case T = 1.5 Tp. 
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FIG. 4: Values of Gp{t,T)/T^ are shown as a function of tT, for T = 3.0 T^. The solid line is the 

result of our calculation made for Gp{T) = Gp[l- 0.17 (T/rj] with Gp = 13.49 GeV-^ while 
the dotted line is obtained when we put Gp{T) = Gp. The data (squares) are taken from Ref. [2] 
for the case r= 1.5 Tc. 
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FIG. 5: Values of Gp{t,T)/T^ are shown for T = 3.0 Tc. Here the solid line is the same as the 
dashed line in Fig. 4 and corresponds to Gp = 0. The dotted line is obtained when we use the 
constant value Gp{T) = Gp = 13.49 GeV^^. 
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solid curve and the dotted curve we see evidence for the temperature dependence we have 
used: Gp{T) ^ Gp[l - 0.17(7/7;) ]■ (We have obtained the value Gp = 13.49 GeV"^ by 
calculating the pion mass at T = in our model.) 

In Fig. 4 we show, as a solid line, our result for Gp(r, T)/T^ at T = 3.0 T^. In this case the 
fit to the data taken from Ref. [2] is poor. Therefore, we have investigated the case in which 
Gp{T) — and show the result as a dashed line in Fig. 4. The fit to the data of Ref. [2 ] is 
improved somewhat. We may suggest that the suppression of Gp{T) at large temperatures 
may be greater than that given by the form we have used, Gp{T) — Gp[l — 0.17{T/Tc) ]. 

In Fig. 5 we again consider the value T = 3.0 Tc. In this case, the sohd line is the result 
obtained when Gp(T) = 0, while the dotted line results when we use a constant value for 
Gp{T) = Gp = 13.49 GeV^^. We may compare our results for the spectral functions to 
those obtained using the MEM procedure and depicted in Ref. [ 2 ] . One essential difference 
is that we do not see the resonant behavior reported there for T — 3.0 at a rather high 
energy of about 2.4 GeV. The peak in ap{uj,T)/uj^ for T = 1.5 Tc is at about 1 GeV in 
Ref. [2], while for T = 1.5 7; our peak in Fig. 1 is at about 0.55 GeV. 

In estimating the energies of the peaks we have used Tc = 150 MeV. However, if we 
use Tc = 270 MeV, which is more appropriate for a lattice calculation without dynamical 
quarks, the resonant structures would appear at still higher energy. 

III. EUCLIDEAN-TIME CORRELATORS OF LORENTZ- VECTOR HADRONIC 
CURRENTS 

We again refer to the Appendix for a review of our procedures for the calculation of the 
spectral functions from which we may obtain the Euclidean-time correlator using Eq. (1.1), 
with av{uJ, T) replacing ap{uj, T). In Fig. 6 we present our values of \mCv{Po- T)/Pq, which 
may be compared to the values of ay (a; , T) /o;^ given in Fig. 1 of Ref. [ 2 ] . In Ref. [ 2 ] we find 
a peak at about 1 GeV for T = 1.5 7; and at about 2.5 GeV for T = 3.07;. On the other 
hand, there is some weak resonant behavior seen for T — 1.5 Tc [dashed hue] in our Fig. 
6 at about 0.75 GeV which may reflect a residual enhancement due to a p-like mode that 
is present after the confinement-deconfinement transition has taken place. At T = 3.0 Tc 
[dot-dashed line], we see no resonance enhancement in our work, in contrast to what is 
obtained by the MEM analysis of Ref. [2]. 
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FIG. 6: Values of ImCviPo,T)/P^, obtained in Ref. [6], are shown for values of T/T^ = 1.2 
[solid line], T/T^ = 1.5 [dashed line], T/Tc = 2.0 [dotted line], T/Tc = 3.0 [dot-dashed line], 
T/Tc = 4.0 [double dot-dashed line], and T/T^ = 5.88 [short dashed line]. Here we use Gv{T) = 
Gv [1 - 0.17 (T/Tc)] with = 0.150 GeV and Gy = 11.46 GeV'^ 

In Fig. 7 we show the data taken from Ref. [ 2 ] as squares. Here, there is hardly any 
difference seen in the data reported for T — 1.5 Tc and T — 3.0 Tc. In Fig. 7, the sohd 
hne represents our results for {?,/ A)Gv{t,T) /T^ at T = 1.5 Tg. We only achieve a fair fit 
to the data, but the fit is decidedly better than that obtained when a constant value of 
Gy{T) = Gy = 11.46 GeV~^ is used [dotted line] (See the appendix for a discussion of the 
factor 3/4 used when making comparison to the data of Ref. [2] in the case of vector current 
correlators.) 

In Fig. 8 we compare our result for (3/4)Gy(T, T)/^^ at T = 3.0 Tc [solid line] with the 
data of Ref. [2 ]. Here the fit is better than that of Fig. 4. Again, the use of a constant value 
of Gv{T) — Gv — 11.46 GeV~^ yields a poor result. We suggest that our analysis tends to 
support our choice of temperature-dependent coupling parameters for the NJL model. 
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FIG. 7: Values of (3/4) Gy(r,T)/T3 are shown for T = 1.5 T^. Here, the soUd hne represents the 

result when wc use Gv{T) = Gy [I - 0.17 (T/Tc)] with Gy = 11.46 GeV"^ The dotted Imc is 
obtained when we use a constant value of Gv{T) = Gy = 11.46 GeV~^. The data (squares) are 
taken from Ref. [2] for the case T = 1.5 Tc. 

IV. DISCUSSION 

It is of interest to obtain further insight into the results shown in Figs. 4 and 5. To that 
end, we show various calculations made for ImCp(Po, T)/Pq in Fig. 9. There, for T — 3.0 Tc, 
the solid line is the result of our model, the dotted curve corresponds to the use of a constant 
value of the couphng parameter Gp{T) — Gp — 13.49 GeV~^, while the dashed hne is the 
result for Gp{T) = 0. A comparison of the solid curve and the dashed curve leads to some 
understanding of the results shown in Fig. 4, while a comparison of the dotted curve and the 
dashed curve leads to further understanding of the results shown in Fig. 5. Similar results 
are given for lm.Gv{PQ,T) / Pq in Fig. 10. For that figure, a comparison of the solid curve 
and the dashed curve gives some insight into the results shown in Fig. 8, where the dotted 
hne corresponds to Gv{T) — Gy — 11.46 GeV~^ and the solid curve represents the results 
of our model. 

We now return to the pseudoscalar case for T = 1.5 Tc. In Fig. 11 we show 
ImCp(Po) 2^)/-Po model [solid line], for Gp{T) = [dashed line], and for the case 
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FIG. 8: Values of (3/4) Gy(r,T)/T3 are shown for T = 3.0 T^. Here, the soUd hne represents the 

result when wc use Gv{T) = Gy [I - 0.17 (T/Tc)] with Gy = 11.46 GeV"^ The dotted Imc is 
obtained when we use a constant value of Gv{T) = Gy = 11.46 GeV~^. The data (squares) are 
taken from Ref. [2] for the case T = 3.0 Tc. 

Gp{T) ^ Gp ^ 13.49 GeV-l We recall that our model, with the temperature-dependent 
couphng parameter, gave rise to a excellent fit to the data, as seen in Fig. 3. It is also of 
interest to present values of ap{u},T) — (l/7r)ImCp(a;, T)/a;^ for T — 1.5 Tc. In Fig. 12, the 
result of our model is shown as a solid line, the dot-dashed hne is for Gp{T) = 0, and the 
dotted line is obtained when Gp{T) = Gp = 13.49 GeV"^ It is seen, that for small values 
of Pq, on the whole, the dotted line lies below the other curves, giving rise to the behavior 
seen in Fig. 3 for tT ~ 0.5. 

It is worth mentioning that we have some additional evidence of the utility of temperature- 
dependent couphng parameters for the NJL model. We note that the confinement- 
deconfinement transition takes place in the range 150 MeV < Tc < 170 MeV for QCD with 
dynamical quarks. We then inspect Fig. 5 of Ref. [ 10 ], where the constituent quark mass of 
the NJL model is presented as a function of temperature for the case of a temperature- 
independent coupling constant. The mass value is 330 MeV at T = and is about 
260 - 300 MeV when 150 MeV < T < 170 MeV. Thus, we do not see the (partial) restora- 
tion of the chiral symmetry that is expected for T ~ Tc. On the other hand, we see in 
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FIG. 9: Values of ImCp{Po,T)/Pg are shown for T = 3.0 T^. The sohd hne is the result of our 

model with tcmpcraturc-dcpcndcnt coupling parameters, the dotted line is obtained in the absenee 
of the temperature dependence {Gp{T) = Gp = 13.49 GeV~^), and the dashed line represents the 
result for Gp{T) = 0. 

Fig. 1 of Ref. [7], where we have used a temperature-dependent coupling parameter, that 
we have m„ = 364 MeV at T = and m„(r) in the range of 50 to 100 MeV for 150 MeV 
< T < 170 MeV. That is much more in accord with the (partial) restoration of chiral 
symmetry when T T^.. If we wish to consider the NJL as a useful low-energy model 
of QCD, it is much easier to discuss the confinement-deconfinement transition if we use 
temperature-dependent coupling parameters. 

For ease of reference,we have calculated the constituent mass of the up quark using the 
equation for the temperature-dependent constituent mass given in Ref. [11]. We have used 
a current quark mass of wP — 5.5 MeV and a momentum cutoff of A = 0.631 GeV. In 
Fig. 13, the dotted curve shows the result obtained with G — 5.691 GeV"^ (in the notation 
of Ref. [11]). From Fig. 13 we see that at T = 0.150 GeV, m„ = 0.318 GeV, while at 
T = 0.170 GeV, m„ = 0.276 GeV. The dashed and solid curves represent the result when 
G{T) = ^[1 - 0.17(T/rc)]. For the solid curve {T^ = 0.170GeV), m„ = 77MeV at T = 
0.170 GeV. For the dashed curve (7; = 0.150 GeV), m„ = 54 MeV at T = 0.150 GeV. 
Again, we see that it is much easier to discuss the (partial) restoration of chiral symmetry at 
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FIG. 10: Values of ImCv{Po,T)/P^ are shown for T = 3.0 T^. The sohd line is the result of our 

model with tcmpcraturc-dcpcndcnt coupling parameters, the dotted line is obtained in the absenee 
of the temperature dependence (Gy(r) = Gy = 11.46 GeV~^), and the dashed line represents the 
result for Gv{T) = 0. 

the confinement-deconfinement transition when we use the temperature-dependent coupling 
parameters of our model. 

It is also of interest to exhibit the role played by the current quark mass. In Fig. 14, the 
dashed curve, which was calculated for Tc — 0.150 GeV, is the same as the dashed curve in 
Fig. 13. In Fig. 14, the solid curve shows the result when m° = 0. Here, we see restoration of 
chiral symmetry at T = 0.136 GeV, when G(T) = G[l - 0.17(T/Tc) ] with = 0.150 GeV. 

We note that the value of Tc = 0.136 GeV is still within the uncertainty of the transition 
temperature for three- flavor QCD. For example, in Ref. [12] the transition temperatures 
are given for two-flavor and three-flavor QCD. In the latter case, T^. — 154 ± 8MeV, with 
a suggested systematic error similar to the statistical error, so that Tc — 154 ± 8 ± 8 MeV 
[12]. (We have used coupling constants determined in our studies of the three-flavor NJL 
model, so that consideration of the transition temperature for that case is appropriate.) 

If we wish to assign a physical interpretation of the parameter Tc in the expression for 
G{T), we may use G{T) = G[ 1 - 0.135(r/Te) ] with Tc = 0.150 GeV. That choice gives rise 
to restoration of chiral symmetry at T = Tc = 0.150 GeV in the NJL model with mP — 0. 
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FIG. 11: Values oilmCp{Po,T)/P^ are shown for T = 1.5 T^. Here, the sohd line is the result of 

our model, the dashed line represent the result for Gp{T) = 0, while the dotted line is obtained 
when we use Gp{T) = Gp = 13.49 GeV-^. [See Fig. 3 for the values of Gp{t, T)/T^ calculated for 
T = l.STc, using the values oilm.Cp{PQ,T) / shown here as the solid and dotted line.] 




FIG. 12: Values of ap{uj,T) = {l/iT)linCp{PQ,T) are shown for T = 1.5 Tc. Here, the solid line 
corresponds to our model, with Gp(T) = Gp[l — 0.17 (T/Tc)], the dot-dashed line is obtained 
when Gp{T) = 0, and the dashed line is for the case Gp{T) = Gp = 13.49 GeV^^. (See Fig. 13.) 
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FIG. 13: We exhibit the values of m„(r) obtained using Eq. (5.38) of Ref. [11], with mP = 
5.50 McV and A = 0.631 GeV. The dotted curve corresponds to the use of a constant value 
G = 5.691 GeV~^, in the notation of Ref. [11]. For the solid and dashed curves we have used 
G{T) = G [ 1 - 0.17(r/rc) ]. For the solid curve we have put = 0.170 GeV, while for the dashed 
curve, we have used Tc = 0.150 GeV in our parametrization of G{T). 

If we maintain the value wP — 0, but use a constant value for G{T) — 5.691 GeV~^, with 
A = 0.631 GeV, we find restoration of chiral symmetry at — 208 MeV. 

The calculations reported in Fig. 13 and 14 may also be made using our Gaussian cutoff, 
exp(— /c ^/a^), with a = 0.605 GcV. For example, we may consider the dashed curve of 
Fig. 14. If we use the Gaussian cutoff and use G = 6.004 GeV~^, instead of G = 5.691 GeV~^, 
we obtain a curve that is very close to the dashed curve of Fig. 14 for T > 0.080 GeV. 
However, the mass at T = is 382 MeV instead of 364 MeV which was the value obtained 
using the sharp cutoff of A = 0.631 GeV and G = 5.691 GeV^^ 

We have performed what are essentially parameter-free calculations of hadronic spectral 
functions and have computed the corresponding Euclidean-timc correlation functions. The 
values of the coupling parameters, Gs and Gy, were fixed in calculations of meson properties 
at T = 0. We have used the temperature dependence, G{T) = ^[1 - 0.17(T/rc)], which 
was introduced in earlier work in which we studied the (mesonic) confinement-deconfinement 
transition [7]. We beheve it is of interest to see that we obtain reasonable values for the 
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FIG. 14: Values of mu{T) are shown. The dashed curve is calculated with mP = 5.50 MeV. Here, 
G{T) = G[l - 0.17 (T/Te)], with G = 5.691 GeV"^ and = 0.150 GeV. The sohd curve is 
calculated with the same value of G{T) and Tc, but with mP = 0. From the solid curve, we see 
that chiral symmetry is restored at T = 0.136 GeV when mP = 0. 

Euclidean correlators. Of more importance, however, is our observation that we find some 
evidence for the temperature dependence of the NJL coupling parameters that we have 
used in other works. By analogy, we expect that the coupling parameters should also be 
density-dependent, and we have introduced such density dependence in earlier work [13]. 
Density dependence of the coupling parameters may be particularly important, given the 
strong interest in diquark condensates and color superconductivity at high baryon density 
[14]. 

APPENDIX 

For ease of reference, we present a discussion of our calculation of hadronic current 
correlators taken from Ref . [ 6 ] . The procedure we adopt is based upon the real-time finite- 
temperature formalism, in which the imaginary part of the polarization function may be 
calculated. Then, the real part of the function is obtained using a dispersion relation. The 
result we need for this work has been already given in the work of Kobes and Semenoff [ 15 ] . 
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(In Ref. [ 15 ] the quark momentum in Fig. 2 is A; and the antiquark momentum is A; — P. 
We will adopt that notation in this section for ease of reference to the results presented in 
Ref. [ 15 ] .) With reference to Eq. (5.4) of Ref. [ 15 ] , we write the imaginary part of the scalar 
polarization function as 

{(1 - ni{k) - n2{k))5{P^ - Ei{k) - E^ik)) 
-{ni{k) - n2(A;))5(P° + E^{k) - E^ik)) 
-Mk) - n,{k))5{P^ - E,{k) + E^ik)) 
-(1 - ni{k) - n2{k))d{P^ + Ei{k) + E2{k))} . 

Here, Ei{k) = [k"^ + m?(T)]^/2. Relative to Eq. (5.4) of Ref. [15], wc have changed the 
sign, removed a factor of g'^ and have included a statistical factor of 2Nc, where the fac- 
tor of 2 arises from the flavor trace. In addition, we have included a Gaussian regulator, 
exp[ —k^/a'^], with a — 0.605 GeV, which is the same as that used in most of our apphca- 
tions of the NJL model in the calculation of meson properties. We also note that 

^iW- eW) + l ' ^^^^ 

and 

= e^^^W + 1 • ^^^^ 

For the calculation of the imaginary part of the polarization function, we may put k'^ — 

m1{T) and {k — P)^ = mKT), since in that calculation the quark and antiquark are on- 
mass-shell. In Eq. (Al) the factor arises from a trace involving Dirac matrices, such 
that 

= -Tr[(|^ + mi)(|^-f + m2)] (A4) 
= 2p2-2(ml + m2)^ (A5) 

where mi and m2 depend upon temperature. In the frame where P = 0, and in the case 
mi — 1712, we have (3s — 2Po^(l — Am^ /Pq). For the scalar case, with mi —m2, we find 

to J,(P^ T) = ^ (l - i^)'" e-«V»'l 1 _ 2n,(k) 1 , (A6) 
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where 



p ^^-m'iT). (A7) 

For pseudoscalar mesons, we replace Ps by 

(5p = -'IV[i75(|^ + mi)i75(|^ - f + ms) ] (A8) 
= - 2(mi - ma)^ , (A9) 

which for mi = m2 is /3p = 2Po in the frame where P — Q. We find, for the tt mesons, 

Im MP\ r) = ^ (l - e-''''-"[ 1 - 2n,(^ ] , (AlO) 

where k"^ — Pq/4: — m^(T), as above. Thus, we see that, relative to the scalar case, the 
phase space factor has an exponent of 1/2 corresponding to a s-wave amplitude. For the 
scalars, the exponent of the phase-space factor is 3/2, as seen in Eq. (A6). 
For a study of vector mesons we consider 

C = 1V[7^(^ + -r + m^)], (All) 

and calculate 

g^'^(3^^ ^A[P'^-ml-ml + Am^m2 ] , (A12) 

which, in the equal- mass case, is equal to APq -(- 8m^(T), when mi = m2 and P = 0. This 
result will be needed when wc calculate the correlator of vector currents in the next section. 
Note that, for the elevated temperatures considered in this work, mu{T) — md{T) is quite 
small, so that AP^ + 9>ml{T) can be approximated by 4Pq, when we consider the vector 
current correlation functions. In that case, we have 

Im Jy (P^ T) ~ ^Im Jp(p2, T) . (A13) 

At this point it is useful to define functions that do not contain that Gaussian regulator: 

~ , , N,P^ ( 4m2(T)\'/' , ,,,, 
Im Jp(P^ T) = ^[\- — ^ j [ 1 - 2ni(/c) ] , (A14) 

and 

Im Jy[P\ T) = [l - — ^ J [ 1 - 2ni(A;) ] , (A15) 
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For the functions defined in Eq. (A14) and (A15) we need to use a twice-subtracted dispersion 
relation to obtain ReJp{P'^,T), or ReJv{P^,T). For example, 

Re Jp(P^ T) = Re Jp(0, T) + [ Re Jp{P^, T) - Re Jp(0, T) ] + (A16) 

P\P' - Pj) f^' , Im Jp{s,T) 

where A^ can be quite large, since the integral over the imaginary part of the polarization 
function is now convergent. We may introduce Jp{P^,T) and Jv{P^,T) as complex func- 
tions, since we now have both the real and imaginary parts of these functions. We note that 
the construction of either Re Jp(P^,T), or He Jv{P'^,T), by means of a dispersion relation 
does not require a subtraction. We use these functions to define the complex functions 
Jp(P^^) and MP^T). 

In order to make use of Eq. (A16), we need to specify Jp(0) and Jp{Pq). We found it 
useful to take P| = -1.0 GeV^ and to put Jp(0) = Jp(0) and Jp{P^) = Jp{Pq)- The 
quantities Jv(0) and Jv{Pq ) are determined in an analogous function. This procedure in 
which we fix the behavior of a function such as ReJy(P^) or ReJv'(P^) is quite analogous 
to the procedure used in Ref. [16]. In that work we made use of dispersion relations to 
construct a continuous vector-isovector current correlation function which had the correct 
perturbative behavior for large P^ ^ — oo and also described that low-energy resonance 
present in the correlator due to the excitation of the p meson. In Ref. [ 16 ] the NJL model 
was shown to provide a quite satisfactory description of the low-energy resonant behavior 
of the vector-isovector correlation function. 

We now consider the calculation of temperature-dependent hadronic current correlation 
functions. The general form of the correlator is a transform of a time-ordered product of 
currents, 

iC{P\ T)^ j dW^- « T{j{x)m) » , (A17) 

where the double bracket is a reminder that we are considering the finite temperature case. 

For the study of pseudoscalar states, we may consider currents of the form jp,i{x) — 
q{x)i^^X^q{x), where, in the case of the tt mesons, i — 1,2 and 3. For the study of scalar- 
isoscalar mesons, we introduce js,i{x) — Qix)^^Q{x), where i = for the flavor-singlet current 
and i = 8 for the flavor-octet current [ 7 ] . 

19 



In the case of the pseudoscalar-isovector mesons, the correlator may be expressed in terms 
of the basic vacuum polarization function of the NJL model, Jp{P'^,T) [11, 17, 18]. Thus, 

1 



Cp{P\T)^Jp{P\T)- 



(A18) 



l-Gp(T)Jp{P\T)' 

where Gp{T) is the coupling constant appropriate for our study of tt mesons. We have 
found Gp{T) = 13.49 GeV~^ by fitting the pion mass in a calculation made at T = 0, with 
TTT'u — TTT'd — 0.364 GcV. The result given in Eq. (A18) is only expected to be useful for small 
P^, since the Gaussian regulator strongly modifies the large behavior. Therefore, we 
suggest that the following form is useful, if we are to consider the larger values of P^. 



cp(p^T) 

p2 



Jp{P\T) 

p2 



1 



Gp{T)Jp{P\T) 



(A19) 



(As usual, we put P = 0.) This form has two important features. At large Pq, 
lmCp{Po,T)/P^ is a constant, since ImJp(Pg^,T) is proportional to P^. Further, the 
denominator of Eq. (A19) goes to 1 for large Pq . On the other hand, at small Pq, the 
denominator is capable of describing resonant enhancement of the correlation function. As 
we will see, the results obtained when Eq. (A19) is used appear quite satisfactory. ( We may 
again refer to Ref. [16], in which a similar approximation is described.) 

For a study of the vector-isovector correlators, we introduce conserved vector currents 
jn,i{x) = q{x)jnXiq{x) with i=l, 2 and 3. In this case we define 

pfj,pi 



4''{P\T)^(g'^^ 



and 



c^(p^T) = u 



p2 



piJ.pv 



p2 



MP ,T) 



Cv{P\T), 



(A20) 



(A21) 



taking into account the fact that the current jij,,i{x) is conserved. We may then use the fact 
that 



MP',T)^-g,^J^^{P',T) 



and 



Im Jy(p2,r) 



P2 + 2ml{T) 



47r 



4ml{T) 



p2 



1/2 



~ -ImJp(P^^). 



(A22) 



[l-2ni(fc)] (A23) 
(A24) 
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(See Eq. (A7) for the specification oi k = \k\.) We then have 

1 



Cy{P\T) = Jy{P\T) 



l-Gv{T)Jv{P\T) 



(A25) 



where we have introduced 



P2 + 2ml{T) 



An 



-\mJp{P\T). 



1-— [l-2ni(A:)] 



(A26) 
(A27) 



In the hterature, ui is used instead of Pq [ 1-3 ]. We may define the spectral functions 

ay (u;, T) = - Im Cy (cj, T) , (A28) 

TT 

and 



ap{u!,T) = — lmCp{u!,T) 



TT 



(A29) 



Since different conventions are used in the hterature [ 1-3 ] , we may use the notation 
ap{uj, T) and ay (cu, T) for the spectral functions given there. We have the following relations: 



ap{uj,T) = ap{u,T) 



(A30) 



and 



ay(a;,T) 3 



CTy{u,T) 



(A31) 



where the factor 3/4 arises because, in Refs. [1-3], there is a division by 4, while we have 
divided by 3, as in Eq. (A22). 



[1] I. Wetzorke, F. Karsch, E. Laermann, P. Petreczky, and S. Stickan, Nucl. Phys. Proc. Suppl. 

106 ,510(2002)-hep-lat/0110132. 
[2] F. Karsch, S. Datta, E. Laermann, P. Petreczky, and S. Stickan, and I. Wetzorke, 

|hep-ph/02Q9028| 



21 



[3] F. Karsch, E.Laermann, P. Petreczky, S. Stickan, and I. Wetzorke, Phys. Lett. B 530, 147 
(2002). 

[4] M. Asakawa, T. Hatsuda and Y. Nakahara, Prog. Part. Nucl. Phys. 46, 459 (2001) 
[5] Y. Nakahara, M. Asakawa and T. Hatsuda, Phys. Rev. D 60, 091503 (1999). 
[6] Bing He, Hu Li, C. M. Shakin, and Qing Sun, |hep-ph/0212345[ 
[7] Hu Li and C. M. Shakin, [hep-ph/ 0209136[ 
[8] J. Letessier and J. Rafelski, |hep^^0301099l 

[9] S. Hamieh, J. Letessier and J. Rafelski, Phys. Rev. C 62, 064901 (2000). 
[10] S. Schmidt, D. Blaschke and Y. L. Kahnovsky, Phys. Rev. C 50, 435 (1994). 
[11] S. P. Klevansky, Rev. Mod. Phys. 64,649 (1992). 
[12] F. Karsch, piep-p h/0106019| 

[13] Hu Li and C. M. Shakin, Phys. Rev. D 66, 074016 (2002). 

[14] For reviews, see K. Rajagopal and F. Wilczek, in At the Frontier of Particle Physics/Handbook 
of QCD, M. Shifman ed. (World Scientific, Singapore 2001); M. Alford, Annu. Rev. Nucl. Part. 
Sci. 51, 131 (2001). 

For a recent application of NJL model in the calculation of the properties of dense matter, 
see D. Blaschke, S. Fredrikesson, H. Grigorian, and A. M. Oztas, |nucl-th/ 03010021 

[15] R. L. Kobes and G. W. Semenoff, Nucl. Phys. B 260, 714 (1985). 

[16] C. M. Shakin, Wei-Dong Sun, and J. Szweda, Ann of Phys. (NY) 241, 37 (1995). 

[17] T. Hatsuda and T. Kunihiro, Phys. Rep. 247, 221 (1994). 

[18] V. Vogl and Weise, Prog. Part. Nucl. Phys. 27, 195 (1991). 



22 



